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Abstract 

The magnetoresistance in the variable-range hopping regime due to Zeeman 
spin-splitting and intra-impurity interactions is calculated analytically and 
shown to be a universal function of fj,H/kT log R. Good agreement with nu- 
merical calculations in one and two dimensions is observed. With the inclusion 
of quantum interference effects, excellent agreement with recent experiments 
is obtained. 
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Magnetic field may have dramatic effects on the resistance of strongly disordered materi- 
als. In this regime, where the resistance is described by Mott law, R = Rq exp{(T /T)*+r}, 
with T ~ l/p£ d , p the density of states, and £ the localization length, the effects of the 
magnetic field can be divided into orbital effects and spin effects. In weak orbital fields, 
of the order of one quantum flux through the hopping area, V-^ 3 £? where L, the hopping 
length, is typically {^/T p) l ^ d+l \ the prefactor of the exponent changes due to quantum 
interference [?,?,?,?], leading usually to negative magnetoresistance (MR). In higher fields, 
of the order of one quantum flux through the area defined by the localization length, £ is 
enhanced [?], leading to an exponential decrease in the resistance. In yet higher fields, the 
impurity wavefunction shrinks, resulting in an exponential enhancement of the resistance 
[?]■ 

As the resistance does not depend explicitly on spin, one may not expect sensitivity of the 
resistance to spin splitting due to magnetic field. However, when intra-impurity interactions 
are taken into account, it was pointed out [?] that the polarization of the electron spins will 
block some of the hopping processes, leading to an exponentially increased resistance. 

It is nontrivial experimentally to separate the contributions of all these mechanisms to 
the MR. This is easier for more disordered samples and at lower temperatures, where the 
magnetic field scales for these effects become more and more separated. More naturally, the 
technological progress in growing strongly disordered thin layers makes it possible to study 
only the spin-effects of a parallel magnetic field [?,?]. It is thus imperative to develop a good 
understanding of the spin effects, not only in order to understand transport in the presence 
of parallel field, but also in order to make it possible to separate the spin effects from the 
orbital effects in more complicated situations. This question became even more important 
with recent suggestions to use the variable-range MR as a sensitive magnetic sensor [?]. 

Kamimura et al. [?], in addition to suggesting the mechanism and deriving the limit of 
high fields, have performed a detailed numerical study of the dependence of the resistance 
on temperature and magnetic field. More recently, Clarke et al. [?] have proposed that the 
resistance may be written, 
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R(H) = R exp{(T /T)^F(y)}, (1) 

with F a universal function of y = pH /T{T /TY^ d+1 . They went on to obtain F(y) in 
limiting situations, and then calculated it numerically as a function of the dimensionless 
parameter y. 

Here I report an exact calculation of the MR due to Zeeman splitting. I find that R{H) 
is indeed given by ([]]), with F(y) obtained exactly in arbitrary dimensions. This function 
agrees very well with the numerical data in one and two dimensions, and with available 
experimental data [?,?]. 

The starting point of this calculation is the mapping of the resistance problem into a 
percolation criterion [?] of an equivalent random resistor network [?], consisting of randomly 
placed sites, of density p, with random energies 6j. The resistance between each pair of sites 
is given by [?] Ry = exp{(|ej| + \cj\ + |e» — £j\)\/2T + 2ry/£}, where all resistances are 
measured in units of Rq. Since the resistances vary exponentially, the overall resistance 
will be dominated by the weakest link, which is the largest resistance, R, such that the 
cluster formed by all resistances (bonds), satisfying R > Rij, percolates. Clearly, all states 
participating in the percolating network (occupied sites) must satisfy R > exp(|ej|/2T). 
Following [?], the percolation criterion employed here is the following [?] - given such an 
occupied site , the probability that a bond is attached to it has to be higher than the critical 
threshold, p c , for the system to percolate, 

Pc = -r^-^ r iOSR &i fde 2 pf d d r 12 Q(R - e (M+M+|e 1 -e 2 |)/2T +2m / e) 

4i lOgR J-2TlogR J J 

2^ 1 r(d/2)d(d + l)(d + 2)' 1 ' 

leading directly to the Mott hopping low. 

In the presence of an intra-impurity interaction U, there are two types of impurities 
participating in the conduction process — those whose energies lie in the vicinity of the 
Fermi energy, Ep (type A), and those whose energy is around Ep — U (type B). For 
U ^> T, the latter are at least singly occupied. Application of a magnetic field, such that 



pH > T, polarizes the singly occupied impurities, thus blocking any process of hopping 
from one singly occupied impurity to another, leading to a lower effective density, and to an 
exponential enhancement of the resistance [?] . For arbitrary magnetic field one can define an 
equivalent resistor network, consisting of sites of two types A and B, with relative densities 
Pa and p B , respectively, such that the resistance between pairs of sites is given by [?,?,?] 
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The percolation condition 
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where pa and p B are the fractions of the sites of types A and B, respectively, on the 
percolating cluster. From ([3]) and (H) one finds that one can write R(H) in the form ([!]) 
with F(x) = {l/[p a p A + PbPb + {PaPB + PBPA)g{x)]} 1/d+1 , with x = 2pH/TlogR, and 
g{x) = 6(1 - x)(l - x) d [l + (d- 2)x/2 + (d - 2)(d - l)a; 2 /12]/(l - x/2). An independent 
equation can be derived for pa and p B , 
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leading to the final expression, F(x) = {2/[l + s(x)]} 1//rf+1 , with s(x) 



(5) 



(Pa — Pb) 2 + 4:g(x) 2 pApB- For small magnetic fields, F(x) ~ 1 + PaPbx, in agreement with 
perturbation theory [?,?]. F(x) saturates at x — 1, or pH = T(T /T m&x{pA, Pb}) 1 +l -> a ^ 
a value of (1/ maxjp^, p B }Y^ d+1 , in agreement with Ref. [?]. 

In Fig.l we compare the exact result to numerical calculation of log R{H)/ log R(0) — 
F(x) in one dimension. In the numerical calculation the resistance of a strongly disordered 
system of length L = 50£ and temperature T = 0.04W 7 , where W is the width of the energy 
distribution, has been calculated using the equivalent random resistor network. There are no 
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free parameters in this comparison. In Fig. 2 we compare the exact result to the numerical 
data in two dimensions reported in [?]. Here the resistance has been calculated using the 
mapping into the percolation problem. Since the parameter x used in that work differs 
from ours and the ratio depends sensitively on the value of To, we allowed a single fitting 
parameter — the x-axis scale. In both Figs. 1 and 2, data has been presented for (a) 
Pa/ Pb = 1, and (b) Pa/ Pb = 1/2. Excellent agreement between the numerical calculations 
and the analytic calculation is observed. 

The real test of the theory is comparison to experimental data. In Fig. 3 we compare 
the analytic result (broken curve) in 2d to the experimental data for the MR of an fr^Os-^ 
layer of thickness d = HOnm in a parallel field [?]. In the inset we compare the analytic 
result in 3d to the experimental data extracted from [?]. In the latter work the magnetic 
field dependence of the the exponent (Eq.[l]) was attributed to an decrease in £ (even though 
theory [?] predicts that £ increases with field). Here the data was replotted in terms of F(x). 
In both plots the two fitting parameters used were pa and the scale of the x-axis. While 
good agreement may be observed in the two figures at high fields, there are clear deviations 
at low fields. This is due to the contribution of the orbital effects, most importantly the 
quantum interference effects. 

The quantum interference effects - the coherent scattering of the hopping electron by all 
other impurities - have been taken into account within the percolation approach in [?]. The 
percolation condition (0) now also involves an interference probability y 2 multiplying the 
exponent inside the function, and an integration over its distribution, which can be derived 
independently [?,?]. The resulting integral equation for the MR has to be solved numerically, 
except for low magnetic fields or deep in the insulating regime. In the latter regime one 
finds that the resistance is multiplied by exp{— < logy 2 >}, where the average is over the 
amplitude distribution. Using direct integration or random matrix arguments, < logy 2 > 
was found to be [?] —7 — log 2 + log(l + J 2 — H/ H^H/ H^) , where 7 is the Euler constant, 
and the field corresponds to one flux quantum through the hopping length. A similar 
procedure can be applied in the presence of spin-split states, as the quantum interference 
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amplitude can be incorporated into the percolation condition (£|) and into Eq.fl^). The 
calculation does not lead to a change in F(x), but leads to a multiplicative factor in Eq. ([[]), 
exp{— A(x) < logy 2 >}, with A(x) = 1 — ApApBxg(x)g' (x) / (d + l)s(x)(l + s(x)). 

Interestingly, this final result suggests that the prefactor resulting from the quantum 
interference is changing even for fields larger than the saturation field H^, due to the change 
in the effective number of impurities participating in the quantum interference. 

The resulting MR is plotted in Fig. 3 (solid line). The three fitting parameters used in 
this plot were pa = 0.76, (To/T) 1 ^ d+1 = 8.2, and H<p = 8.1 Tesla (no fit was used to scale the 
x-axis). The latter two values should be compared to the approximate experimental values, 
(Tq/T) 1 / 3 ~ 10.6, and (po/dL ~ 9.1 Tesla, respectively. An almost perfect agreement with 
the experimental data is observed, especially in view of the fact that in the experimental 
paper [?] a 5-parameter fit was used, and even then the agreement, it was claimed, "by itself, 
was not trivial" , within the physical constraints on the parameters used. Similar procedure 
can be successfully applied to fit the data of Ref. [?] (inset of Fig. 3). 

For higher temperatures or perpendicular magnetic fields, the last step in the calcula- 
tion - the mapping of the percolation problem into a log-averaging procedure, is no longer 
applicable [?]. The formalism is still valid, but one needs to solve the full integral equation 
discussed above. 

In this work we have reported an exact calculation of the magnetoresistance due to 
spin effects in the variable-range-hopping regime. The calculation compares very well with 
available numerical and experimental data. The mechanism responsible for that effect — the 
blocking of hopping from a singly occupied impurity to another singly occupied impurity due 
to spin polarization suggests that spin-orbit scattering may play an intriguing role in such 
systems. As the temperature is lowered the hopping length gets larger and may eventually 
become larger than the spin-orbit scattering length. When this happens the electron may flip 
its spin upon hopping from one impurity to another, thus making it possible to hop from one 
singly occupied impurity to another. Thus it is predicted that spin-orbit scattering may lead 
to nontrivial temperature effect and even to a decrease in the resistance as the temperature 
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is lowered. As such "turning-on" of spin-orbit scattering with lowering of the temperature 
has been demonstrated for the magnetoresistance in the weakly localized regime [?], an 
experimental investigation of this prediction should not be too difficult. We hope that this 
work will motivate such investigations. 
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Figure Captions 

1. Comparison of the numerical calculation of log it! (H)/ log R (0) in one dimension to 
the analytic F(x), with no free parameters. 

2. Comparison of the analytic F(x) to the numerical data extracted from Ref. [?] in two 
dimensions. The only free parameter is the x-axis scale. 

3. Comparison of the analytic result with (solid line) and without (broken line) including 
interference effects to the experimental data of [?] and (in the inset) of [?]. The fitting 
parameters are discussed in the text. 
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